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We present the results of the spectroscopic observations of HDE 226868, the optical 

1 ^ . counterpart to the black hole X-ray binary Cyg X-l, from 2001 to 2006. We analyze 
the variabilities of the two components in the complex Ha line: one P-Cygni shaped 

component which follows the motion of the supergiant and another emission component 
moving with an antiphase orbital motion relative to the supergiant, which is attributed 
to a focused-stellar wind. The results of KOREL disentangling of our spectra indicate 
that the focused stellar wind is responsible for the major part of the variability of the Ha 
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2 ■ emission line. The emission of the supergiant component had a small difference between 

the low/hard and high/soft states, while the focused wind component became strong 
in the low/hard state and weak in the high/soft state. The wind is nearly undisturbed 
by the X-ray photoionization during the low/hard state. However, during the high/soft 
state, the X-rays from the compact object could decelerate the line-driven wind and 
' result in a high mass accretion rate, due to the effect of the X-ray photoionization. The 

CN ■ X-ray illuminating could also change the temperature profile of the stellar wind and 

increase its temperature, and thus decrease the Ha emissivity of the wind, which could 
explain the Ha variabilities of Cyg X-l during different X-ray states. 
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Subject headings: Stars: Binaries: Spectroscopic, Stars: Early- Type, stars: individual 
(HDE 226868, Cygnus X-l), Stars: Winds, Outflows, X-Rays: Binaries 



INTRODUCTION 



Cygnus X-l was first discovered during a rocket flight observation in 1964 (IBowver et al- 
and its optical counterp art was identified with the supergiant star HDE 226868 (jBolton 



!rp3 
3). 



1965) 



1972 



Webster & Mur dinlll972T ). A 5.6d orbital period was found in the optical spectroscopic observations 



(jGies &: Boltonlll982l ) and the dynamical determination of the binary components provided evidence 



for the existence of a black hole in Cyg X-l. As the black hole X-ray binary, which was first found 
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in our Galaxy, Cyg X-l has been extensively studied from radio t o 7-r ays during recent decades. 
Even though the mass function is precisely known (see lGies et al.l (|2003l )). the masses of the donor 
star and the black hole are not so well constrained due to the poorly know n inclination and the 
evolutionary status of the supergiant. Using the inclination of i = 35°, iHerrero et al.l ()1995l ) 



estimated the masses of the supergiant and the black hole at 17.8 M Q and 10.1 M , respectively. 



HDE 226868 was classified as an 09.7 lab supergiant star (|Walbornlll973l ). which shows Ha 
and He II A4686 in emission. Although the supergian t nearly fills the Roche lobe, the accre tion is 
mainly via the strong stellar wi nd from the donor star (I Gies k Bolton 1986a ; Gies et al. 2003 ) . The 
variabilities of the optical lines ( Gies k Bolton 1986b ; Ninkov et al. 19871 ) on the spectrum of HDE 
226868 indicate that the distribution of the stellar wind deviates fro m a spherical geometry and 
that an enhanced wind flow exists ("focused stellar wind," suggested by lFriend k Castorl (|1982l )) in 
the direction of the compact object. The focused stellar wind is also revealed by the High-Ener gy 
Transmission Grating Spectrometer aboard the Chandra X-Ray Observatory (jMiller et al.ll2005l ). 



The most recent ep hemer is of the 5.6d orb ital period has been given by lLaSala et al.l ([1998), 
Brocksopp et al.l (|1999bl ) , and IGies et al.l (|2003l ) according to their optical spectroscopic observa- 
tions. In addition, the 5 .6d orbital period is also found in the UBVJHK photometry, X-ray and 
radio data of Cyg X-l (jBrocksopp et al.1 Il999al ). Cyg X-l also shows the superorbital modula- 
tion, on a time-scale mu ch longer than the or bital period. A 29 4d period was first reported by 
Priedhorsky et al.1 (|1983l ) in X-rays an d then bvlKemp et al.l (119831) in the optical. An other ~ 150d 



period was found by different authors (jBrocksopp et al.lll999al : lLachowicz et al.ll2006l ) in X-ray and 
radio data. This ~ 150 d period may be caused by the precession of the accretion disk around the 



compact object (jWiiers k Pringldll999l ). 



A relativistic radio jet with a velocity larger than ~ 0.6c was detected bv lStirling et al.l (|200ll ) 
in their radio observations of Cyg X-l, and therefore Cyg X-l joins the group of the Galactic 
microquasar. Microquasar is an X-ray binary with a p air of relativistic radio j ets, w hich is similar 
to the radio jets found in the active galactic nuc lei fsee Mirabel k Rodriguez! (119991) for a review) . 
The jet remnant of Cyg X-l is resolved in radio (|Gallo et alJhoOfJ ) and optical (IRussell et al.ll2007l ) 
observations. 

Many observational characteristics of the canonical stellar black hole candidate Cyg X-l are 
considered evidence that a black hole exisits, similar to other X-ray binary systems. These observa- 
tional features include the ultra-soft spectra, the high-energy power-law tail above 20 keV, the spec- 
tral/temporal transition from th e high / soft state to the low /hard state and the X-ray millisecond 



variab ility (see iTanaka k Lewinl ()1995l ): iMcClintock k Remillardl (j2006l ): iRemillard k McClintock 



(J2006J)). Most of the time Cyg X-l is in a low/hard state and in some years it can transit from 
the low/hard state to the high/soft state, w hich can continue from several weeks to several months 



(jZhang et al. 



1997 



Brocksopp et al.lll999al ). Occasionally, Cyg X-l enters an inter mediate state 



when it fails to make a transition from the low/hard state to the high/soft state (jBelloni et al 



1996 



Malzac et al.ll2006l ). The radio and X-ray emission has an anticorrelation and when the bi- 
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nary enters the high/soft state, the radio emission will be decreased Most researchers 

believe the transition sta t e is c aused by the physical changes of the gas around the black hole (see 
Remillard &: McClintockl (|2006l ) and references therein). The physical changes in the accretion disk 
are initially related to the mass loss from the supergiant companion. However, the X-ray radiation 
of the disk influences the ionization and temperature of the wind and thus also its radiation and 
dynamics. 

Using the method of tomographic separation^ ( Bagnuolo h Gies 1991 ). Sowers et al. ( 1998 ) 
showed that the Ha profiles of Cyg X-l observed in seasons 1985 and 1986 can be well fitted as a 
superposition of a P-Cyg profile corresponding to the (approximately spherical) stellar wind of the 
supergiant and a wide emission peak radiated by the focused stellar wind, which moves in a slightly 
shifted anti-phase with respect t o the orbital Doppler shift of the former component. Using the 



method of Fourier disentangling (jHadravalll995l . 



1997 



2004 ) we have proved the same behavior in 



the spectra obtained at Ondfejov Observator y in summer 2 003 shortly before, during and shortly 
after one high/soft-state episode of Cyg X-l (jHadraval 120071 ). 



In this article, we present our optical spectroscopic observations of Cyg X-l from 2001 t o 2006. 
The preliminary results have been introduced in a previous article (jYan. Liu &: Hang||2005l ). Here 
we will make a further analysis of the Ha line profiles. Our present data are confined to about one 
week in each season, meaning they are concentrated into only slightly more than one orbital period. 
This enables us to compare the solutions of disentangling in different X-ray states, minimizing the 
possible influence of long-term variations in the structure of the focused stellar wind. 

In the following section, the properties of our data are described. The method and results of 
the disentangling are given in Section 3. Next, the correspondence of our observational results with 
the standard model of this classical X-ray binary is studied in Section 4. Finally, the conclusions 
are summarized in Section 5. 



2. OBSERVATIONS 

We obtained the spectra of HDE 226868 with the 2.16m telescope at Xinglong Station of 
National Astronomical Observatories, China (NAOC), from 2001 to 2006. The optical spectroscopy 
with an intermediate resolution of 1.22 A pixel -1 was made with a CCD grating spectrograph at 
the Cassegrain focus of the telescope. We took the red spectra covering from 5500 to 6700 A and 
blue spectra covering from 4300 to 5500 A at different times. Sometimes low-resolution spectra 
(covering from 4300 to 6700 A) were also obtained. The journal of our observations is summarized 
in Table [H including observational date, UT Middle, exposure time, Julian date, wavelength range, 



1 This method of decomposition of observed spectra into unknown spectra of two components should be dis- 
tinguished from the methods of Doppler tomography in which the phase-locked line-profile variations are fitted by 
projections of smooth distribution of delta-function profiles in velocity space corotating with the orbital motion. 



- 4 - 



and spectral resolution. Orbital phase (0) is also given in Table [JJ and the ephemeris of the inferior 
conjunction of the companion star is adopted from Gies et al. (2003), 

2, 451, 730.449 + 5.599829£ . 

All spectroscopic data were reduced with the IRAfU package. They were bias-subtracted and flat- 
field corrected, and had cosmic rays removed. Helium-argon spectra were taken in order to obtain 
the pixel-wavelength relations. To improve this relation, we also used the diffuse interstellar bands 
(DIBs) 6614 and 6379 A observed in the spectra. 

The higher resolution red spectra obtained from 2001 to 2006 are shown in Figure [TJ The corre- 
sponding observational dates and orbital phases calculated according to the above given ephemeris 
are written on the left and right sides of each spectrum, respectively. The spike on the left part of 
the spectrum on 2002 October 26 may be caused by a cosmic hot point and the dips on the spectra 
of 2004 September 21 are caused by bad pixels on the CCD. Most of the Ha lines show a double- 
peaked profile with a central absorption. Single-peaked Ha lines are observed in our spectra of 
2003 October 14 (0=0.67), 2005 October 24 (0=0.004) and 2006 September 29 (0=0.72). Obvious 
P-Cygni Ha lines are observed in our spectra in some phases only in 2004 (Figure [JJc)). 

For a comparison, we have also used the spectra obtained with the 700mm camera of the Coude 
spectrograph at the 2.05m telescope of the Ondfejov observatory (the Astronomical Institute of 
the Czech Academy of Sciences). These spectra covering t he region 6260-67 60 A with a resolution 



of approximately 0.25 Apixel 1 are included in a study bv lGies et al.l ( 20081 ). to which we refer for 
details. 



3. ANALYSIS AND RESULTS 



3.1. The Equivalent Width Evolution of Ha and the X-ray Activity 



The equivalent width (EW) of the complex Ha line (emission and absorption) has been mea- 
sured selecting a continuum point on each side of the line and integrating the flux relative to the 
straight line between the two points using the procedures available in IRAF. The measurements 
were repeated five times for each spectrum and the error estimated from the distribution of values 
obtained. The typical error for Ha measurements is within 10%. This error arises due to the 
subjective selection of the continuum. The results of Ha EWs are listed in Table [TJ 

The top panel of Figure [2] shows the variability of the Ha EW from 2001 to 2006 as a function 
of time. In addit ion t o the data obta i ned in our observational program, the combined data sets of 
Gies et al. (j2003l ) and Tarasov et al. ( 2003 ) are also included in the figure. For a comparison, the 



2 IRAF is distributed by NOAO, which is operated by the Association of Universities for Research in Astronomy, 
Inc., under cooperation with the National Science Foundation. 
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middle panel of Figure [2] gives the RXTE/ASM one-day averaged counter rate in the 1.5-12 keV 
band and the hardness ratio (HR1) of the soft X-ray radiation, (3-5 keV)/( 1.5-3 keV) is plotted in 
the bottom panel. The arrows in this panel correspond to the starting date of each observational 
run. 

Figure [2] shows that our 2001 and 2004 observational runs were done in a high/soft state 
while 2003, 2005, and 2006 observational runs in a low/hard state. The observations in 2002 were 
during a transitional state from high/soft to low/hard. The EW of Ha is relatively low in the 
hi gh /soft state and stro ng in the low/hard st ate. This phenomenon has been discussed in detail 
bv lGies et al.l (|2003l ) and lTarasov et al.1 (|2003l ). While the Ha in our spectra of 2001 nearly lost its 
emission signature, the emission level of Ha in the 2006 observations is the strongest among our 
six observational runs. Because the season 2001 is poorly covered by observations, we chose the 
observational runs 2004 and 2006 to represent the high/soft and low/hard states, respectively, in 
our study of line-profile variability. 



3.2. The Profile Variability of Ha Emission Line 



The top panel of Figure [3] shows the Ha profiles during our 2006 observational run. The 
exposures are depicted in ascending order according to the orbital phase. The observational date 
and orbital phase are marked on the left and right sides of each spectrum, respectively. We have 
recalibrated the wavelength scale of each spectrum according to the position of the DIB 6614 A. 
It can be seen from this figure that, for most of the time, the Ha of Cyg X-l has a double- 
peaked profile, one peak formed in the supergiant and the oth er in the focused stellar wind between 
the system components (powers et al.lll998l : iGies et al.ll2003h. As in many other supergiants, its 
intrinsic Ha emission is due to the powerful stellar wind (IPuls et al.lll996l ) and it forms the red- 
shifted emission wing of the P-Cygni line profile that follows the orbital motion of the star. In some 
phases (e.g., 0.720, which is close to one extreme of the radial velocities), this emission merges with 
the emission of the focused wind moving almost in an anti-phase and they form a single bright 
peak. In some other phases, the emission of the focused wind can fill the absorption part of the 
P-Cygni component formed in the supergiant. Thus we can detect a characteristic P-Cygni Ha 
line in the spectrum of Cyg X-l only in some phases of the high/soft state, when the focused- wind 
emission is weak. The bottom panel of Figure [3] shows a gray-scale map of the Ha profiles in 2006. 



Following ISowers et al.l (|1998l ). lGies et al.l (|2003l ) used the tomographic separation algorithm to 
decompose the Ha line profile into the two components for chosen combinations of the radial velocity 
semiamplitude K em and phase shift <p of the focused wind (orbital parameters of the supergiant 
being fixed from nonemission lines). Minimizing the residuals in the two-parameter space, they 
found that the focused wind component has a maximum radial velocity K em = 218±30kms _1 near 
the orbital phase </>o = 0.79=L0.04. The radial velocity curves of the two components are plotted in 
Figure HI The positions of the focused-wind component of the Ha line based on this solution are 
also marked by ticks in Figure [TJ 
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3.3. Disentangling the Ha Line Profile 



The methods of disentangling fit the observed spectra as a superposition of several components 
with simultaneously optimized orbital parameters (and/or some other free parameters). Usually, 
the emission-line objects are variable on a time scale shorter than the orbital period. This creates 
problems for disentangling, tomographic separation, Doppler imaging or any other similar method 
that requires observations spanning a time interval at least of this order (unless the variability is 
properly involved in the model of component spectra). However, there is always a chance that such 
methods can reveal a mean behavior of the object treating the rapid variability as a noise. 



The results by ISowers et al.l (|1998l ) and iGies et al.l (|2003l ) indicate that this assumption can 
work well in the case of Cyg X-l. To decompose the spectra of the supergiant and the focused stellar 
wind, we thus use the KOREL code (|Hadravall2004l ) for Fourier disentangling (|Hadravalll995l ). which 
enables us t o take into acc ount and to resolve instantaneous changes in stren gth of lines of each 
component (IHadraval 1 199 71 ) . In the preliminary results for the Ondfejov data (jHadraval 120071 ) the 
He I line at 6678 A has been disentangled into a weak telluric contribution and a pure absorption 
line of the supergiant, which confirms within the observational errors the orbital parameter s, foun d 
by classical me asurements ( e.g., K\ = 71.9kms _1 compared to 75.6 ± .7 by IGies et al.1 ([2003), 
or 73.0 ± .7 by IGies et al.l (|2008l )). The Ha line has been disentangled into three components 
corresponding to the supergiant, the focused wind, and the telluric water vapor lines, which are 
relatively strong in some Ondfejov spectra. Because, unlike the orbital velocities of the supergiant 
and the black hole, the mean velocity of the focused stellar wind is not perpendicular to the line 
joining the components of the binary, the option of KOREL to disentangle up to five components in a 
hierarchical structure has been used to identify the focused wind with a component of a second close 
pair corotating with the supergiant-black hole pair with identical period but with a shifted phase. 
The results showed the P-Cyg profile moving with the supergiant and a broad emission peak of the 
focused wind. The orbital parameters of the P-Cyg profile converged in the disentangling to values 
consistent with the He-line solution (e.g., K\ = 71.3 km s -1 ). The variations of its strength were 
relatively small (of the order of 0.1); moreover the EW of this component is also small, because the 
red emission wing nearly compensates the blue absorption wing. For the broad emission component 
we found a s e miam plitude of K em = 60.8 km s -1 , whic h agrees bet t er wi th the value 68kms _1 by 
Sowers et al.l (|1998l ) rather than the value obtained by IGies et al.l (|2003l ); compare their figures 5 
and 6, respctively. The broad minima of the spectra residuals in the parameter space are due to 
the large width and the variability of the focused-wind emission component, and the position of 
the deepest point may be influenced by long-term changes in the circumstellar matter as well as 
by the random sampling by the observations. The line-strength factor derived by disentangling 
of the Ondfejov spectra for the focused wind is significantly higher for this component (reaching 
a value around +1 at the initial and final low/hard Ha emitting states). Because the absolute 
value of EW of the focused wind (which is negative) is higher, the variability of this component is 
responsible for the major part of the enhancement of the Ha emission in the low/hard state. In 
the disentangling of the Ondfejov data, a small part of the Ha emission appeared in the telluric 
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spectrum, which was left arbitrary to test the behavior of the solution. Because the annual changes 
of the heliocentric radial velocity corrections are appreciably smaller compared with the amplitudes 
of both the supergiant and the focused-wind disentangled components, this part of the emission 
belongs to circumstellar matter which is in the mean at rest with respect to the center of mass of 
the binary system. 

In disentangling the NAOC spectra we struggled with two instrumental obstacles. One was 
the smaller spectral resolution and the other was the unreliability of the dispersion curve, probably 
caused by insufficient stiffness of the Cassegrain spectrograph. Due to the former, no telluric lines 
can be seen in these spectra. Fortunately, it seems that they are not as strong as they are in some 
Ondfejov spectra, in which they can be distinguished even after a smoothing. Consequently, they 
do not need to be disentangled to clean the stellar spectra, but at the s a me ti me they cannot be 



used to check or improve the wavelength scale, as described by lHadraval (|2006l ). We thus tried to 
use the above-mentioned DIBs 6614 and 6379 A for this purpose. We measured their positions 
in the Ondfejov spectra first to find their mean central wavelengths in Cyg X-l and then in the 
NAOC spectra to get the same values by a linear transformation of the wavelength scale in each of 
these spectra. For the measurement we used a sin gle-component d isentangling by KOREL in the 



option of free radial velocities in each exposure (e.g. . lHadraval (|2004l )) . This procedure improved the 
wavelength scale to some extent, yet a solution of the He I line 6678 A made in order to check the 
reliability of the improved wavelength scale revealed errors in some exposures up to about 60 km s _1 
and similarly the disentangling of Ha did not provide satisfactory results. The insufficient precision 
of the results may be caused by the weakness of the DIB 6379 A, which was still measurable in the 
Ondfejov data, but too shallow and wide in the lower-resolution data from the NAOC. 

We thus chose an alternative way to disentangle the Ha profiles in spite of unreliable wavelength 
scale. First we disentangled both seasons 2004 and 2006 together with the Ondfejov data for the 
three components (the supergiant, the focused wind and the telluric lines) in the option of free 
velocities. The velocities and line strengths were kept fixed for the Ondfejov exposures from their 
previous standard solution. This solution was thus used as a template to which radial velocities and 
line strengths of each NAOC exposure were adjusted. The line strengths were also converged by 
the simplex method (instead of direct least squares which would also change the Ondfejov values) 
for the first two components in these exposures and prescribed to a large negative value for the 
telluric component to diminish it in the NAOC exposures. (The difference in spectral resolution 
of these two data sets proved to be unessential in these wide profiles.) The results of this solution 
were used as the initial approximation for two component solutions (without the telluric lines) 
of free radial velocities and line strengths (now by the direct calculation) independently for both 
seasons 2004 and 2006 without the other data. These solutions thus cannot yield correct radial 
velocity curves for each component and proper wavelengths of the disentangled spectra, but they 
give correct differences between the radial velocities of the two components, their line profiles, and 
the line strengths. 



The disentangled Ha profiles and the fit of the input spectra for our 2004 and 2006 observations 
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are plotted by the standard graphical KOREL output in Figure [5] at the left and right panels, 
respectively. These spectra span from 6540 A to about 6596 A; however, the shift of each profile 
up to several A is uncertain due to the above-explained problems. It can be seen (at the upper 7 
and 11 curves, respectively.) that the agreement of the fit with the observation is quite good (some 
differences appear for the fifth and sixth exposures from the top in the right panel only, i.e. for 
2006 September 28). 

The bottom two curves show the disentangled P-Cyg profile of the supergiant (the higher 
of these two curves) and the mean disentangled emission profile of the focused stellar wind (the 
very bottom curve) in each season. The tops of the emission wings of the P-Cyg profiles are 
1.06 and 1.09 of the level of continuum in 2004 and 2006, resp. This indicates a 50% increase 
above the continuum in the low/hard state in 2006; however, this result deserves a confirmation 
by more extensive observations, because the small change representing 3% of the continuum may 
be influenced only by uncertainties of the continuum intensity (e.g., due to differences in phase 
sampling). In any case, this change in the P-Cyg component is negligible in comparison with the 
overall enhancement of Ha emission in the low/hard state. The absorption wing is shallower in 
2004 than in 2006 with the deepest point at levels 0.93 in 2004 and 0.83 in 2006. It also can be 
seen that the decrease of the emission intensity toward the higher velocities is somewhat slower in 
2006 and similarly the depth of absorption is more pronounced for higher velocities in the 2006 
low/hard state. The EW 0.54 A of the absorption part of the P-Cyg profile in 2004 is canceled 
(with precision almost 10 -3 A) by EW -0.54 A of the emission part of the profile. In 2006, the 
EW of P-Cyg absorption is +1.910 A and emission -1.916 A. Note that these EWs refer to the 
disentangled mean seasonal profile, for which the noise is decreased and continuum is fitted to an 
extended spectral region. Consequently, their errors are smaller than those given in Table [T] for 
individual exposures (which are in the mean ±0.03 A and ±0.05 A in 2004 and 2006, respectively). 
The continua of disentangled components may suffer from some complementary distortions induced 
by errors in rectification of the input spectra. The number of significant digits given here is to show 
that emission and absorption EWs of the P-Cyg component cancel each other out in both states 
within the precision of our data. 

More remarkable difference is found between the disentangled profiles of the emission of the 
focused wind: its maximum is 0.056 in 2004, and 0.265 in 2006 in the units of the continuum level 
of the supergiant (a possible continuum of the focused wind cannot be disentangled, but it must 
be negligible compared with the continuum of the supergiant). Similarly, the mean EWs of the 
focused-wind emission are -0.33 A in 2004 and -1.78 A in 2006. 

The line-strength factors, like the EWs, seem to show some phase-locked variations, which 
could reveal the geometry of the system (i.e., the distribution of the focused stellar wind in the 
space between the components) and eventually also the anisotropy of the deeper layers of the 
stellar wind in the upper atmosphere of the supergiant. However, the amplitude of these changes 
is comparable to the observational scatter and the phase coverage of our data is not yet sufficient, 
so we postpone this problem to our next study. 



- 9 - 



Our results generally confirm that it is the contribution of the focused wind, and not the P-Cyg 
profile formed in the root of the supergiant's wind, that is responsible for the major part of the 
Ha emission enhancement in the low/hard state compared with the high/soft state. The physics 
of this process will be discussed quantitatively in the following section. 



DISCUSSIONS 



4.1. The X-ray Excited Wind Model 

The wind from isolated O and B stars is accelerated to a speed of approximately three times the 
escape velocity from the surface of the star by the force arising from the absorption and scattering 
of the photospheric continuum radia tion in the ultraviolet resonance lines of abundant ions in the 
wind ([Castor. Abbott fc Kleinlll975l ). The radiation driving force produces a velocity profile of 



v (r) 



,(1.0 -R*/r) 



(1) 



in a steady state, where v(r) is the wind velocity at a distance of r from the center of the star, 
oo is the terminal velocity of the stellar wi nd, i?* is the radius of the supergiant, and (3^0.8 



([Friend fc Abbottlll98q ; |Pauldrach et al.lll986l ). The gravitational effect of the compact companion 
in a massive X-ray binary system induces a stream of enhanced wind (focused stellar wind) in the 
line from the supergiant to the compact ob ject. The leve l of th e densi ty enhancemen t in th e focused 
wind is less than a factor of 2 as derived bv lHaberl et al.l ([ 19891 ). while iBlondin et al.l (jl99ll ) thought 
that the density in the focused wind was 20-30 times the ambient wind density. The presence of 
the focused wind in the massive X-ray binary could greatly enhance the mass accretion rate of 
the compact object in the system. Meanwhile, the X-ray radiation from the compa ct object can 



strongly in fluence the dynamics of the wind via X-ray heating and photoionization (jBlondin et al 



1990 



1991 



The acceleration of the wind can be inhibited by the X-ray photoionization which can enhance 
the degree of ionization in the stellar wind. Thus, the deceleration of the wind due to the pho- 
toionization will greatly enhance the mass accretion of the compact object, which is a sensitive 



function of the wind velocity law, M a 



~(XV 



and lead back to a higher X-ray luminosity. The 



temperature and ionizatio n state of the stellar wind depend only on the "ionization parameter" 
(|Kallman fc McCravl[l982h . 

£ = L x /n p r 2 x , (2) 

where L x is the X-ray luminosity of the compact object, n p is the proton number density of the 
wind, and r x is the distance to the X-ray source. When £>10 2 ergs cm s _1 , the stellar wind 
will be strongly affected by X-r ay ionization and w hen £>10 3 ergs cm s _1 , X-ray heating will 
affect the dynamics of the wind (jBlondin et al.lll990l ). To estimate the ionization parameter, let 
us suppose the wind of the supergiant to be spherically symmetric with a constant mass-loss rate, 
M=4irr 2 nm p n p (r)v(r) , where fi is the mean atomic weight and m p is the mass of the proton. Then 
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the ionization pa rameter in terms of the position and the stellar wind parameter can be derived 
(|Sako et alJll99gh . 

C(r,r x ) = 4.3 x 1Q 2 i^M ^^ _ ^ (3) 



M. 



7 



' x 



where {L x )^ is the X-ray luminosity in a unit of 10 36 ergs s , (foo)s is the terminal velocity in a 
unit of 10 8 cm s _1 , and M_7 is the mass loss of the supergiant in a unit of 10 Mq yr _1 . 



Using the transformation formulae of IZdziarski et al.l (j2002l ). we can convert the RXTE/ASM 
one-day averaged c ounts in three band s (1.5-3, 3-5, and 5-12 keV) to the energy units. Adopting a 
distance of 2.5 kpc (jNinkov et al.lll987l ) to the source, we can get the X-ray luminosity of Cyg X-l 
in differe nt X-ray states. T h e mas s-loss rate of the 09.7 supergiant is calculated by the formulae 
given by iHowarth Sz Prinial (|1989l ). M=2.OxlO _6 M yr~ x . All parameters of Cyg X-l are listed 
in Table The ionization parameters during our 2001, 2003, 2004 and 2006 observational runs 
are calculated according to Equation ([3]) for two distances r x =21R & and 23i?0 from the black hole 
and they are listed in Table El We also plot the contours of constant £ (£=10 2 ) for Cyg X-l in 
different X-ray states in Figure El where the coordinate origin is in the center of the supergiant, 
the thick circle is the surface of the supergiant, and the intersection of the two black dashed lines 
is the boundary of the supergiant's Roche lobe. 



4.2. The Mass Accretion Rate in Different X-Ray States 

The X-ray activity of the compact object is closely related to the mass accretion rate onto 
it. Since the supergiant component of Cyg X-l has a negligible change in the Ha line between 
the low/hard and high/soft X-ray states, as discussed in the previous section, the mass-loss rate 
from the supergiant does not have an obvious change between these two types of X-ray states. 
One possible factor causing the enhancement of the mass accretion rate during the high/soft state 
may be the X-rays from the compact object. Figure [6] indicates that the ionization parameter £ 
has a value of about 100 ergs cms -1 at the position of r x =21R & (in the line from the compact 
object to the supergiant) during the high/soft state, while it only has a value of 11 ergs cm s -1 
during the low/hard state. Consequently, the focused wind may be photoionized by the X-ray 
during the high/soft state, which leads to the decrease of the radiative pressure exerted in lines 
by the supergiant's own radiation, and hence also the velocity of the gas will be decreased. The 
slower wind velocity will greatly enhance the mass accretion of the compact object and a high 
X-ray luminosity will be observed during the high/soft state. During the low/hard state, the X-ray 
radiation is not strong enough to affect the focused wind, hence Cyg X-l has a relatively low-mass 
accretion rate. 

A related question we should discuss is whi ch kin d of m echanism causes Cyg X-l to transit 



from the low /hard state t o the high /soft state. iDond (|2002l ) suggested that the disk instability 



mechanism (DIM; lLasotal (|200ll )) triggered the X-ray outburst and then the X-ray irradiation 



photoionized the hydrogen in the wind of the primary and enhanced the mass accretion onto the 
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compact object. The occurrence of the X-ray outburst does not have any periodicity and it may 
happen at any instant in the course of time. The duration of the outburst may be connected with 
the interaction between the X-rays from the compact object and the wind of the supergiant. In 
general, it can keep several weeks or even several months (see Figure [2]). 



4.3. The X-Ray Irradiation Effects on the Ha Emission 



X-ray heating coul d affect the dyn amics of the wind from the supergiant. According to the 
plane-parallel model bv lWu et al.l (|200ll ). a temperature-inversion layer will be formed in the stellar 
wind under the X-ray illumination. The temperature profile of the temperature-inversion layer is 
determined by the soft and hard X-ray irradiation. A strong temperature-inversion layer can be 
formed provided that the incident X-rays are soft and their angle of incidence is approaching grazing 
incidence. The more penetrating hard X-rays tend to heat up the deeper layers of the stellar wind. 

During the high/soft state, the increased soft X-ray heating will increase the temperature 
of the ambient wind gas and then decrease the Ha emissivity, which is proportional to T~ L2 
([Richards Ratlin 11998). The X-ray heating is more effective in the position near the X-ray 
source so the focused wind component in the Cyg X-l is more easily affected by the X-ray radiation 
from the compact object. Thus, the focused-wind component in the Ha line becomes weak during 
the X-ray high/soft state. During the low/hard state, the temperature profile of the focused wind 
is undisturbed by the X-ray irradiation and thus it has a strong emission in the Ha line. Compared 
with the focused wind, the supergiant is farther from the X-ray source and its spherical wind is 
rarely influenced by the X-ray heating. However, when the X-ray luminosity is very high, the wind 
near the supergiant also can be affected (see Figure [6]). For the high/soft state during our 2004 
observational run the outer parts of the supergiant stellar wind were more ionized and thus unable 
to absorb and re-emit in Ha. The terminal velocity of the wind also was smaller. For the low/hard 
state during our 2006 observational run, the soft X-rays from the compact object have almost 
no influence on the supergiant stellar wind. This could explain the small differences between the 
shapes of disentangled profiles in the two states in Figure [5j When the soft X-ray flux was very 
strong, the wind near the supergiant could also be influenced by the X-ray photoionization. Thus 
the system might lose its emission feature both in the supergiant and the focused-wind components. 
This scenario is consistent with the weak emission Ha line in our 2001 spectra (on the bottom of 
Figure (Ha)). 

During our 2004 observation, Cyg X-l was also in a high/soft state, but the X-ray intensity 
was not so strong as that in our 2001 observations. Therefore, the focused-wind component in the 
Ha line was at a lower emission level, while the supergiant component had a small change relative 
to the case in the low/hard state. Then we cou l d obs erve a n obvious P-Cygni H a profile during our 
2004 observations. The spectra of iGies et al.1 (|2003l ) and iTarasov et al.l (|2003l ) also indicate that 
the P-Cygni Ha profile often appeared in the X-ray high/soft state. During the low/hard state, we 
rarely observed the P-Cygni structure in the complex Ha line of Cyg X-l, because the focused-wind 
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Ha component is strong enough to fill the absorption part of the P-Cygni line formed in the wind 
of the supergiant. 



4.4. The X-ray Flare During Our 2006 Observations 

The Ha EWs in our 2006 observational run are stronger than those in other runs (see Figure [2]). 
Figure Old) shows that the focused-wind component in our 2006 Ha spectra becomes very strong 
while the supergiant component is still in its normal emission level. Since the supergiant component 
did not have a significant change during our 2006 observations, it is inappropriate to explain the 
enhanced emission in the focused-wind Ha component using the increased mass-loss rate from the 
supergiant. Figure [2] indicates that there is a small X-ray flare around the 2006 observational 
run. This X-ray flare was also d etected by Swift /BA T in the 15-20 keV high-energy band. This 



phenomenon has been reported by I Albert et al.l (|2007l ). Figure [7] shows the Ha EW evolution during 



our 2006 observations and the light curves of the flare detected by RXTE/ASM and Swift/BAT, 
respectively. The dashed lines correspond to the beginning and the ending times of our 2006 
observations. 

This flare differs from the X-ray outburst during the high/soft state. It has a relatively low 
X-ray flux and a h ard X-ray spectrum. The TeV emission also was detected during the flare by 



Albert et al.l (|2007l ) and they suggested that the flare was caused by the interaction between the jet 
and the stellar wind. Our observations were carried out on the decline phase of this flare. The weak 
soft X-ray emission during the flare is insufficient to form a prominent temperature-inversion layer 
in the focused-wind of Cyg X-l and therefore a strong focused wind Ha component was observed. 
However, the hard X-ray component during the flare could heat up the deeper layer of the focused 
wind and an extra Ha emission could be detected during our spectroscopic observations. This can 
explain the strongest Ha emission during our 2006 observations. 



Another similar strong Ha line with an EW of -2.315 A was also observed by iTarasov et al 



(|2003l ) on MJD 5,0941.5449, when a small X-ray flare (low/hard state) was also detected by 



RXTE/ASM in the 1.5-12 keV band (see the inset in the top panel of Figured]). 



5. CONCLUSIONS 

We present and analyze our optical spectroscopic observations of Cyg X-l from 2001 to 2006. 
Combined with the RXTE/ASM X-ray observations, we make the following findings about this 
classical galactic black hole X-ray binary: 



1. We confirm that the Ha line shows two components: a P-Cygni profile moving with the radial 
velocity curve of the supergiant and a focused stellar wind component that moves with an 
approximately anti-phase orbital motion relative to the supergiant. The superposition of the 
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two components forms the complex Ha profiles. 

2. The results of KOREL disentangling the Ha spectra during our 2004 and 2006 observations 
indicate that the focused-stellar wind is responsible for the major part of the Ha variability 
between different X-ray states. The focused wind component becomes strong during the X- 
ray low/hard state and weak during the high/soft state. The photoionization and heating of 
the X-ray photons from the compact object may affect the ionization state and dynamics of 
the wind from the supergiant. During the high/soft state, the X-ray photoionization could 
decelerate the gas in the focused stellar wind and result in an increasing mass accretion rate. 
During the low/hard state, the X-ray is not strong enough to influence the wind ionization 
state and the compact object has a low-mass accretion rate. The X-ray illumination can form 
a temperature-inversion layer in the stellar wind. During the high/soft state, the soft X-rays 
acting onto the focused stellar wind could increase its temperature greatly and thus decrease 
the Ha emissivity. This could explain the variability of the focused-wind component in Ha 
during different X-ray states. The strong soft X-ray emission during the high/soft state could 
also ionize the outer parts of the supergiant winds and render it unable to absorb and re-emit 
in Ha. This scenario is consistent with the small differences between the shapes of disentan- 
gled P-Cygni components in the 2004 and 2006 spectra. During our 2001 observations, the 
wind near the supergiant was also affected by the strong X-ray emission and an extremely 
weak Ha line was observed. 

3. The Ha lines in our 2006 observations are very strong. We interpret this as the result of 
the low irradiation of the focused wind by the X-ray photons. The weak incident soft X-rays 
during the flare could not disturb the temperature profile of the focused stellar wind, while 
the hard X-rays could heat up the deep layer of the wind and an extra Ha emission could be 
observed. 

J.Z.Y. is grateful to Min Fang for his help in plotting Figure [3] and Figure [6j This research is 
partially supported by the National Natural Science Foundation of China under grants 10433030 
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Fig. 1. — Spectra including Ha and He I A6678 lines during six years' observations. The observa- 
tional dates and orbital phases are given in the left and right sides of each spectrum, respectively: 
(a) Ha spectra in 2001 and 2002, (b) Ha spectra in 2003 and 2005, (c) Ha spectra in 2004, and 
(d) Ha spectra in 2006. When the Ha shows a double-peaked profile, the focused-wind component 
is marked by a tick in the figure. 
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JD-2,450,000 



Fig. 2. — (a) Ha EWs during six- year observations (filled circles). Data from iGies et al.l (J2003) 
(open circles) and iTarasov et al.l (j2003l ) (open triangles) are shown, (b) The one-day-averaged 
RXTE/ASM count rates of Cyg X-l in 1.5-12 keV. (c) The hardness ratio in the soft X-ray band 
of RXTE/ASM, (3-5 keV)/(1.5-3 keV). The arrows on the bottom of the panel correspond to the 
starting time of each observational run. 
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Fig. 3. — Top: Ha profiles during the 2006 observations, arranged in ascending order according to 
the orbital phase. The observational dates and orbital phases are given on each side of the Ha, 
respectively. All spectra have had the continuum level normalized and offset vertically to allow 
direct comparison. Bottom: the gray-scale map of the Ha spectra in 2006. The y-axis corresponds 
to the orbital phase and the gray intensity is scaled between 0.986 (white) and 1.17 (black). 
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Fig. 4. — Radial velocity curves of th e supergian t (soli d line) and focused-wind (dotted line) com- 
ponents in the Ha line adopted from iGies et al.1 (120031 ) . The two components have a comparative 
velocity around orbital phases 0=0.03 and 0.53. 
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Fig. 5. — Disentangled Ha lines observed in high/soft state in the 2004 (left panel) and low/hard 
state in 2006 (right panel). The upper 7 lines (left panel) or 11 lines (right panel) indicate the 
input spectra ordered in running time from the top. These profiles are superimposed with their 
reconstruction as the sum of disentangled components, which are shown by the two bottom curves. 
The very bottom lines correspond to the emission of the focused stellar wind, which is obviously 
higher in 2006 than in 2004. The lines second from the bottom display the P-Cyg profiles of the 
supergiant. The maxima of their emission wings are practically equal in both states. 
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Fig. 6. — Contours of £=10 2 for Cyg X-l in different X-ray states. The coordinate origin is at 
the center of the supergiant and the thick circle represents the surface of the supergiant. The 
intersection of the two black dashed lines is the boundary of the supergiant's Roche lobe. The 
solid, dashed, dot-dashed, and dot-dot-dashed lines are the contour lines when the X-ray luminosity 
(L x ) 36 equals 14.28 (2001), 11.47 (2004), 1.46 (2003), and 2.04 (2006), respectively. 
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Fig. 7. — The top panel shows the Ha EWs during the 2006 observations. An X-ray flare was 
detected by Swift/BAT (15-50 keV) (middle) and RXTE/ ASM (1.5-12 keV) (bottom) around our 
observations. The dashed lines correspond to the starting and ending times of the 2006 observations. 
The inset in the top panel is another X-ray flare detected by RXTE/ASM (1.5-12 keV) around MJD 
50,941.5449 (the po sition of the arrow), when a strong Ha emission (-2.315 A) was observed by 
Tarasov et all ifcooah . 



Table 1. Summary of the spectroscopic observations of HDE 226868. 
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1000 


2453275. 


,0635 


5550- 


■6750 


1 


.22 





.832 





,28±0. 


,01 


20040926 


13 


:49: 


:02 


1000 


2453275. 


.0757 


5550- 


■6750 


1 


.22 





.835 


0, 


25±0, 


,01 


20051021 


13 


:30 


:33 


900 


2453665. 


0629 


5550- 


•6750 


1 


.22 





,477 


1 


,27±0. 


.04 


20051023 


11 


:31: 


;50 


100 


2453666. 


.9804 


4300- 


■6700 


2 


.44 





.820 


1, 


53±0, 


,06 


20051023 


11 


:34: 


:05 


100 


2453666. 


.9820 


4300- 


•6700 


2 


.44 





.820 


1. 


.54±0. 


.01 



Table 1 — Continued 



Date UT Exposure Julian Wavelength Spectral Orbital Ha 

Middle Time Date Range Resolution Phase EW 

(hh:mm:ss) (s) (A) A pixel -1 (-A) 



20051024 


12: 


;17: 


40 


900 


2453668. 


.0119 


5550- 


-6750 


1 


,22 


0. 


,004 


1 


.04±0. 


,02 


20051024 


12: 


:32: 


:50 


900 


2453668. 


,0228 


5550- 


-6750 


1 


.22 


0. 


.006 


1 


.08±0. 


,02 


zUUolUz / 


1Z. 


.00, 


OU 


con 


z40OD 1 1. 


Uoyo 


4oUU- 


-oouu 


1 
1 


.zz 


n 
U. 


KA K 
.040 








20060926 


14: 


:05: 


:43 


1200 


2454005. 


,0873 


5550- 


-6750 


1 


.02 


0. 


.198 


1 


,75±0. 


.07 


20060926 


14: 


:27: 


48 


1200 


2454005. 


,1023 


5550- 


-6750 


1 


.02 


0, 


.200 


1 


,79±0. 


.06 


20060927 


13: 


:26: 


:20 


1200 


2454006. 


,0600 


5550- 


-6750 


1 


.02 


0, 


.372 


1 


,71±0. 


,02 


20060927 


14: 


:05: 


:07 


1200 


2454006. 


,0869 


5550- 


-6750 


1 


.02 


0, 


.376 


1 


.74±0. 


.07 


20060928 


12: 


:24: 


:02 


1200 


2454007. 


,0167 


5550- 


-6750 


1 


.02 


0, 


.542 


1 


,63±0. 


,09 


20060928 


13: 


:25: 


:46 


1200 


2454007. 


,0596 


5550- 


-6750 


1 


.02 


0, 


,550 


1 


.76±0. 


,08 


20060929 


12: 


:19: 


:54 


1200 


2454008. 


,0138 


5550- 


-6750 


1 


.02 


0, 


.720 


1 


,96±0. 


,03 


20060930 


13: 


:28: 


1 1 


1200 


2454009. 


,0616 


3850- 


-5050 


1 


.02 


0, 


,908 








20061001 


12: 


:04: 


:59 


600 


2454010. 


,0035 


5550- 


-6750 


1 


.02 


0, 


.076 


1 


.69±0. 


,02 


20061001 


15: 


:09: 


:09 


600 


2454010. 


,1314 


5550- 


-6750 


1 


.02 


0, 


099 


1 


.56±0. 


,03 


20061002 


12: 


:09: 


:59 


600 


2454011. 


,0069 


5550- 


-6750 


1 


.02 


0. 


.255 


1 


.95±0. 


,04 


20061002 


14: 


:35: 


:00 


600 


2454011. 


,1076 


5550- 


-6750 


1 


.02 


0. 


.273 


1 


.96±0. 


,03 



The ephemeris is adopted from Gies et al. (2003). 



Table 2: Parameters of Cyg X-l 



Parameter 


Value 


Reference 


M* 


24±5 M Q 


Iorio (2007) 


M x 


8.7±0.8 M 


ShaDoshnikov & Titarchuk f2007) 


q=M x /M, 


0.36±0.05 


Gies et al. (2003^1 


a 


42±9 R Q 


Iorio (2007) 


i 


48.0±6.8° 


Iorio (2007) 


Roche lobe size 


r m =21±6 R G 


Iorio (2007) 




1700 km s -1 






40,000K 




d 


2.5 kpc 


Ninkov et al. ( 19871 


L pt/ Lq 


5.59xl0 5 


Ziolkowski f2005) 


M w 


2.OxlO" 6 M yr -1 


Howarth k Prinia C1989) 


R* 


18 R Q 





-26- 



Table 3: Ionization parameters in different X-ray states 



Observational 


X-ray- 


X-ray luminosity 


£(^=21^0) 


C(r x =23R Q ) 


Run 


State 


(10 36 ergs s^ 1 ) 


(ergs cm s _1 ) 


(ergs cm s^ 1 ) 


2001 


High/Soft 


14.28 


110.04 


40.89 


2003 


Low/Hard 


1.46 


11.25 


4.18 


2004 


High/Soft 


11.47 


88.38 


32.84 


2006 


Low/Hard 


2.04 


15.72 


5.84 



